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Abstract 

We investigate homogeneous and isotropic oscillating cosmologies with multiple fluid 
components. Transfer of energy between these fluids is included in order to model the 
effects of non-equilibrium behavior on closed universes. We find exact solutions which 
display a range of new behaviors for the expansion scale factor. Detailed examples are 
studied for the exchange of energy from dust or scalar field into radiation. We show 
that, contrary to expectation, it is unlikely that such models can offer a physically viable 
solution to the flatness problem. 

1 Introduction 

In this paper we consider a wide range of homogeneous and isotropic cosmological models con- 
taining two fluids which are able to exchange energies and manifest non-equilibrium behavior. 
Recently [1] , we studied a class of such examples in the dynamical context of a flat Friedmann 
universe. A simple ansatz was used to model energy transfer between two fluids, which would be 
separately conserved perfect fluids in the absence of the energy transfer. By means of this sim- 
ple formulation we were able to find a single master equation whose range of solutions allowed 
us to study the cases of massive particles decaying into radiation, particle-antiparticle anni- 
hilation into radiation, the decay of classical vacuum energy, and the formation and Hawking 
evaporation of black holes. In this paper, we add some further ingredients: the effects of spatial 
curvature are included in the cosmological dynamics and a wider range of energy exchanges 
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are considered. The most interesting situation allowed by the introduction of spatial curvature 
is the possibility of an oscillating closed universe. Following Tolman's [2] identification of the 
important role that the second law of thermodynamics could play in the long term evolution of 
a cyclic universe in which entropy production takes place, there have been further detailed stud- 
ies of how thermodynamics leads to differences in the cycle-to-cycle evolution. Landsberg and 
Park [1] studied the approach to flatness created by the growing cycles; Barrow and Dabrowski 
[5] carried out a detailed study of the evolution of anisotropies and black holes, and also showed 
that if a positive cosmological constant is present then Tolman's growing oscillations will even- 
tually be ended - replaced by unending expansion towards the de Sitter space-time. However, 
these studies did not attempt to model the action of energy exchanges between matter and 
radiation; they simply injected entropy (in the form of radiation) into the universe at each 
moment of bounce. A more realistic modelling of energy exchange was introduced by Barrow, 
Kimberly, and Magueijo [6], who included energy transfer between scalar fields and radiation 
in some cosmological models arising in general relativity, Brans-Dicke theory, and theories in 
which variations in the fine structure 'constant' are driven by appropriately coupled scalar fields 
[TJ. This revealed the interesting feature that oscillating universes, in which time variations in 
the gravitational 'constant' and the fine structure 'constant' can occur, saw those variations 
continuing in an almost monotonic fashion from cycle to cycle. Unlike the expansion scale fac- 
tor, these varying 'constants' did not oscillate. A further variant was considered by the present 
authors, in ref [8], where they considered energy exchange between radiation and the Brans- 
Dicke scalar field so that energy was drained from the gravitational 'constant' as the universe 
expanded. Here, we extend these investigations and enlarge the gallery of possible cosmological 
evolutions of the scale factor that can occur in oscillating closed universes in general relativity. 



2 Cosmological Models 



Consider a spatially homogeneous and isotropic Friedmann universe, with expansion scale factor 
a(t), containing two fluids with equations of state 

V = (7 - 
Pi = ( r - !)Pi> 

where the 7 and T are constants, p and p\ are the fluid pressures and p and p\ are their densities. 
The Friedmann equation is then given by 

H 2 = P + Pl -- 2 , (1) 
a 1 

where 87rG/3 = 1, H = a/a is the Hubble expansion rate of the universe, and t is comoving 
proper time. We write the evolution equations for these two fluids as 

p + m 1P = 8 , (2) 

p 1 + 3HTp 1 = -s (3) 
where s is a function parameterising any exchange of energy between them. 
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In general, s could be any function of H, p, p 1 ,t and a. In a previous work pQ, we considered 
the case where s was prescribed by 

8 = -aHp + pHpi, (4) 

with a and /3 both constant, and k — 0. This case can be used to model the two-way exchange of 
energy between fluids in a number of physically interesting situations involving particle decays, 
vacuum decay, particle-antiparticle annihilation, or black hole formation and evaporation, and 
is particularly appealing because it can be solved completely in terms of simple functions. Here, 
we will consider more general formulations of this problem, particularly those that arise when 
there is non-zero spatial curvature (A; ^ 0) and with different dependence of the energy exchange 
parameter, s, on the other physical variables. 

The inclusion of spatial curvature allows particular cases of special interest, such as oscillat- 
ing universes with entropy production, to be investigated. In these models a universe of positive 
spatial curvature repeatedly expands to a maximum and recollapses towards a big crunch. In 
the absence of non-adiabatic processes, these cyclic universes will all have the same maximum 
of expansion and total lifetime. However, it was first argued by Tolman [2] that an increase of 
entropy at the moment of each crunch-to-bang transition should create a growth in the scale 
of successive maxima of expansion. In this way, the universe may approach spatial flatness at 
asymptotically late times, allowing a potential solution of the flatness problem. These growing 
oscillations will continue until such time as any positive cosmological constant, or quintessence 
field that violates the strong energy condition, dominates the expansion, after which the oscil- 
lations will first shown in ref [5]. If the entropy production occurs in relatively small 
steps then the final state should also be one in which the energy density in the quintessence 
field is just slightly greater than that in the cold dark matter fields. 



3 Cosmological Thermodynamics 



The exchange of energy and momentum between different components of a universe containing 
multiple fluids requires us to reconsider the thermodynamics of the universe, especially the 
effects of the second law of thermodynamics. We write the fundamental law of thermodynamics 

as 

TdS = d{ P V) + pdV = d[(p + p)V]-Vdp (5) 



where S (not to be confused with s) is the entropy of the fluid and V 
volume of the universe. The integrability condition 

d 2 S d 2 S 



is the comoving 



dTdV dVdT 



then gives the relation 



dp 



(P + P) 
T 



dT 



which can be substituted into equations (j2J) and (jSJ to give 



d_ 

dt 



~ {p + P)V ~ 
T 



Vs 



—— and 



dS = d 



~ {p+p)v ~ 

T 
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We can then read off the expression for the density p fluid 



S 



sV 



(6) 



and similarly, for pi, 

S, = -f • (7, 

J 1 

It can be seen directly from these expressions that allowing a transfer of energy between different 
components of the universe (s ^ 0) results in the entropy of those components no longer staying 
constant. In the limit that this interaction disappears, s — > 0, the rate of change of entropy also 
disappears and entropy is again conserved. These expressions for the rate of change in entropy 
of the two fluids enable us to assess the implications of the second law of thermodynamics. 
According to the second law: 



This equation tells us that for energy to be transferred from p to pi (i.e. s < 0) we must have 
T ^ Ti and that, conversely, for energy to be transferred from pi to p (s > 0) we must have 
T\ ^ T. In other words: energy is always transferred from the hotter to the cooler fluid. 

4 Exact Solutions 

There are a number of exact solutions that can be found to the system of equations ([1]), ([2]) 
and 0. In this section we will present methods of integrating these equations, the analysis of 
which will follow in the subsequent sections of the paper. 

4.1 Solutions with p = p(a) 

Taking the form of s given by (J3J) and making the definitions 



S + S 



sV 




A = a + /3 + 3( 7 + r), 
B = ^(aT + /3 7 + 3 7 r), 



the equations ([2]) and ([3]) can be put in the simple form 



— + A— + 25 = 

P P 



where primes denote differentiation with respect to the new time coordinate, 



rj ee In a. 
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The solution to this equation is then given by 

where c\ and c 2 are constants. Substitution into ([3]) then gives pi oc p and so ([1]) can be 
rewritten as 

H 2 = c 3 a^- A -^^) + c 4 ai(- A +^™) - A- 

a 2 

It can be seen that this equation reduces to the solution found in ref pQ in the limit k — > 0, and 
that the inclusion of spatial curvature does not add any new features to the form of a that are 
not already present in the standard evolution without interaction. Specifically, it can be seen 
that in this case an oscillating universe with k > will have the same maximum of expansion 
in each of its cycles. 



4.2 Solutions with p = p(a, t) 

The Friedmann equation ([I]) shows that, for any functions p = p(a) and p\ = pi(a), the 
evolution of the universe shares the same property as the example above. Oscillating k > 
universes have constant amplitude. In order to find new behaviors, we will consider a more 
general situation where we parameterise p = p(a, t) and p\ = pi(a, t). This will allow new exact 
solutions to be found with two interacting fluids, and non-zero spatial curvature. Integrating 
equations (j2J) and (j3J), we obtain 

f sa 3l dt m — f sa 3T dt mi 
P = ^— = -§r and pi = ^ = 8 

which defines the quantities m and mi. In the absence of energy exchange (s = 0), these two 
quantities are constant and we have the same form for p and pi as in the usual case. As soon 
as energy exchange is allowed, m and m x become non-constant, corresponding to energy being 
exchanged between the two fluids. 



4.2.1 A universe containing radiation and dust or massless scalar field 

For the cases of universes containing radiation and dust, or radiation and a scalar field, the field 
equations ([1]), (T5]) and ([3]) can be transformed into the equation of motion for a forced harmonic 
oscillator. Such equations can be solved exactly using standard techniques. These cases are of 
particular interest for modelling an oscillating universe as they correspond to physical significant 
situations. It has often been hypothesised that an oscillating universe containing these fluids 
should allow for energy to be transferred, and entropy to be increased, at the moment of crunch- 
to-bang. The case of exchange between radiation and a scalar field is of particular interest as 
the kinetic parts of scalar fields may dominate the earliest stages of the universe, and non- 
equilibrium behavior can arise through slow or fast decays, as expected at the end of inflation 
or other phase transitions involving these fields. Negative-energy scalar fields are often used in 
oscillating universe models to produce a non-zero minimum of expansion, and a detailed study 
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of their effects, especially on any 'constant' of Nature which may vary in time was given in ref 
®. 

Differentiation of the Friedmann equation (jTj), and some manipulation and substitution 
using (|2J) and (J3J) allows us to write the evolution equation for H(t) as 

* + !r* = ?(r- 7 )„ + ( i-!r£ 

Following the prescription given in ref [9] to turn non-dissipative Friedmann universes into 
simple harmonic oscillators in conformal time, we make a transformation to the conformal time 
coordinate adr = dt and define the new variable b = a 37-2 , to obtain the equation of motion of 
a forced harmonic oscillator: 

6" + (2-3 7 ) 2 A;6= i(2-3 7 ) 2 m, (9) 

where primes here denote differentiation with respect to r. In deriving this equation we have 
used the definition of m given by (JBJ), and assumed the relation 3r + 2 = 67: This includes 
the important cases of radiation and dust (r = 4/3 and 7 = 1) and radiation and scalar field 
(7 = 4/3 and T = 2), which this section focuses upon. The well known solution to equation 
(JSj), for general m and k 7^ 0, is given by 

6(r) = ci sin (V(2 - 37) 2 A;r) + c 2 cos (V (2 - 37) 2 fcr) 

^ (2 2 ~ 37)2 sin (V(2-3 T ) 2 fer) J mcos (V(2 - 3^n) d n 

- v/(2 2 ~ 37)2 cos(v / (2-37) 2 A:r) ^ m sin (V(2 - 3 7 ) 2 £:t 1 ) dn, 

where ci and C2 are constants of integration. By specifying a particular energy-exchange func- 
tion, s, the value of m can be found from equation (jBJ), and the integrals above can then be 
performed. Methods of obtaining various forms of m(t) are outlined in the appendix. 

If we take, for example, the forms of s to be 

sa = SoP a 1 an d sb = ■s a~ ( ' 1+37 ' ) , 

where So is a constant, then the corresponding functions m are found in the appendix to take 
the particularly simple forms 

rriA = m e s ° T and tub — rn + s r, 

where mo is a constant of integration. The solutions, for k 7^ 0, can then be written as 



+ 



oa = C\ sm 

and 



c\ sm 



2((2 - 3 7 ) 2 A; + si 

V¥^Kr-f Q ))^^±^ (11) 

where c±, c\, r and f are constants, and subscripts A and B denote solutions corresponding to 
sa and sb, respectively. The behavior of these solutions will be investigated in the next section. 
The case of no energy being exchanged can be read off from equation fill I) as the sub-class of 
solutions with s = 0. 
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4.2.2 A universe containing radiation and non-zero vacuum energy 



A second mechanism for obtaining an oscillating universe is a negative vacuum energy. Due 
to the lack of dissipation of vacuum energy as the universe expands, this fluid always comes 
to dominate at late-times (if the universe survives that long, and in the absence of phan- 
toms). Negative vacuum energy then plays the role that was previously played by positive 
spatial curvature, and causes the universe to collapse at late-times as anti-de Sitter behavior is 
approached. 

We will now proceed to find solutions for a universe containing radiation and non-zero 
vacuum energy (p = —p). Again, the Friedmann equations can be cast into the form of the 
equation of motion for a harmonic oscillator, but now the forcing term will be constant and the 
analogue of the mass will be non-constant. 

Taking T = 4/3, 7 = and defining the new variable a = a 2 , allows the equations ([I]), (j2J) 
and (j3J) to be reduced to the form 

4A 

a - —a = -2k, (12) 

where use has been made of (|HJ), and A = 3m = 3p has been defined in analogy to the usual 
notation of the cosmological constant. (The reader should bear in mind that A is not a constant 
here, as energy is being exchanged between it and the radiation fluid). This equation can be 
solved once s has been specified, and a solution obtained for A. 

If we now take, for example, the energy-exchange parameters 

sc = s A and sc = s 

then, from the appendix, we obtain the functions 

A c = A e 3s °* and A D = A + 3s t, 

where so and Ao are constants. The equation (fT2l) can now be solved for these s and m, to give 

2k f 2k f 

ac = C3/0 (x) + c 4 A (x) - —I Q (x) / K (xi) dtx + — K (x) / J (x x ) dt x (13) 

3s J 3s J 

and 



a D = c 3 A i {y)+c i B i {y) + ks Q ^I^A l {y) j Bi{y x ) dh - ks ^j ^ B i(v) J A{yi)dh (14) 

where C3, C3, C4 and C4 are constants, Jo and K are Bessel functions, A4 and Bi are Airy 
functions and x and y are defined by 



X=W W e3S ° </2 ^ V= V2^I (A ° + 3S0t) ' 



The form of a given in ({TBI corresponds to A > 0; making the substitutions K — > nJ and 
I — > Y gives the solution for A < 0. For k = the unsightly terms involving integrals of 
Bessel and Airy functions in (fT3l and (IT4"|) vanish. 
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In the absence of any energy exchange, A is constant, and the solution to equation (1121) is 

. ( /-4A . .1 3k 
« = c 5 sm|^^(t-t 1 )j + ^ 

where C5 and ti are constants. In the next section we will proceed to investigate the physical 
behavior of these solutions. 

5 Behavior of the Cosmological Solutions 

Having found exact solutions for oscillating universes with energy exchange, in the previous 
section, we will now investigate their behavior. The oscillatory nature of these solutions is 
created either by positive spatial curvature, or a negative vacuum energy, either of which will 
cause a maximum of expansion, after which the universe collapses. By introducing energy 
exchange between the fluids in these models, we allow for the possibility of different cycles 
having different expansion maxima. 

5.1 A closed universe containing radiation and scalar field 

From equations ({TO]) and (flT|) above, we see that the scale-factor for a universe containing 
radiation (7 = 4/3) and scalar field (T = 2) with positive spatial curvature (k = 1) can be 
written as 

4 = Cl sin{2(T-7- )} + e so(T - ri) (15) 

when s = Sopa" 1 , and as 

a| = Cisin{2(r - f )} + ^s (t - fi) (16) 

when s = sqci~ 5 . Here we have absorbed some of the constants of ( 1101) and (TlTl) into the new 
constants T\ and f\. 

The two functions fTTS]) and ffTBI) above have the undesirable feature of allowing a 2 to be 
both positive and negative, for any set of the constants ci, r and t\. However, only two of 
these three will be fixed by specifying the energies of the two fluids at the moment of collapse; 
the third is a free constant which can be used to match together different cycles, at the moment 
of crunch-to-bang. In this way the scale-factor can be made real and positive semi-definite, 
throughout its evolution. Figures [1] and [2] show the evolution of these universes after such a 
matching. 

It can be seen from figures [JJ and [2] that, as time progresses, the expansion maximum of 
each cycle increases. This is because we chose so > to construct these plots, corresponding to 
energy being transferred from the scalar field to radiation. This causes successive maxima of 
expansion to increase in magnitude, due to the extra radiation. This is in accord with Tolman's 
original model [2]. However, as well as the maxima increasing from cycle to cycle, we also see 
that after a certain number of cycles there exist real non-zero minima of expansion, instead of 



S 




T 



Figure 1: The evolution of the scale factor in a closed universe with radiation and scalar field 
exchanging energy, as prescribed in equation (TB) when < s oc pa -1 . 




Figure 2: The evolution of the scale factor in a closed universe with radiation and scalar field 
exchanging energy, as prescribed in equation /[To}) when < s oc a~ 5 . 

collapse to a singularity. After this point the evolution of the universe is non-singular. It can 
clearly be see that these minima of expansion also increase in magnitude as r increases. This 
is an effect that has not been recognised in earlier discussions of oscillating universes subject 
to the second law of thermodynamics. We interpret this behavior as being due to the scalar 
field transferring so much energy to the radiation field that its own energy density becomes 
negative, allowing a non-zero minimum of expansion. Negative energy scalar fields are often 
used in this way to create oscillating universe models without singularities; however, we can 
include a cut-off to the energy transfer that would prevent this from occurring. 

In addition to this common behavior, significant differences can also be seen to occur between 
figures [T] and [2] at both early and late-times. Such differences in behavior clearly illustrate that 
the evolution of these universes is sensitive to the particular form of energy exchange that has 
been chosen. 

At late-times, figure [1] can be seen to exhibit a runaway behavior. As the energy density of 
the radiation in this universe increases, so does the rate of energy transfer (as it is proportional 
to p). This results in the endless transfer of energy, at ever increasing rates, from the unbounded 
negative energy scalar to the radiation. In fact, after a sufficient number of cycles the runaway 
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transfer of energy eventually becomes so great that it overwhelms the oscillatory nature of the 
solutions, after which these universes expand eternally. The late-time behavior of figure [2] also 
allows an endless transfer of energy from the negative energy scalar to radiation. However, now 
the energy transfer is no longer increased by the ever-increasing energy density in radiation, p, 
and so does not display the runaway behavior seen in figure [TJ The increase of each expansion 
maximum and minimum happens at a steadier rate, and at no point does the energy exchange 
become rapid enough to overwhelm the oscillatory nature of the solutions. 

At early-times, these solutions exhibit cycles which are separated by catastrophic collapse 
to a singularity. As time is run backwards, the maximum of each cycle decreases, as the energy 
of the radiation decreases. Figure [1] shows approach to a regime where each cycle is of the 
same amplitude, in the limit, as p — > and hence sa — > 0. Expansion is then due to the scalar 
field, until it is inevitably halted by the positive spatial curvature. By contrast, in figure [2] we 
see that the expansion maximum and the total duration of each cycle continues to decrease, 
as time is run backwards. This behavior occurs because the energy density of radiation is no 
longer bounded from below, and continues to decrease until it becomes negative. After this 
point, the radiation contributes to the onset of collapse, increasingly so as its energy density 
becomes more negative. Of course, negative energy radiation may not be particularly realistic, 
and we can include a cutoff in the energy transfer function to prevent this from occurring. 
Such a cutoff will eventually result in cycles of constant amplitude, when energy is no longer 
being transferred. It might also be avoided by picking forms for s which are more physically 
motivated, and may have a more complicated evolution in time, rather than the simple choices 
we have made for illustrative purposes. 

The case of energy transfer from radiation to the scalar field can be pictured using figures 
[Hand [2] by reversing the direction of time. 

5.2 A closed universe containing radiation and dust 

From equations ffTUl) and ffTTj) we now see that a closed universe containing radiation (r = 4/3) 
and pressureless dust (7 = 1) has scale-factor 

a A = ci sin(r - f ) + e so(r " fl) (17) 

when s = s pa _1 , and 

a B = cisin(r - f ) + -s q (t - h) (18) 

when s = soa~ 4 . Again, we have absorbed some of the constants of ( fTOl) and (llip into the new 
constants f% and f\. 

The behavior of the scale-factor is qualitatively the same as that shown in figures [T] and 
121 (with the axis label a 2 being replaced by a) . The main difference here is that s > (the 
direction of time going from left to right in the figures above) now corresponds to energy 
being transferred from radiation to dust. This results in the subsequent maxima and minima of 
expansion increasing as the energy density of the dust increases. For the more realistic situation 
of energy going from dust to radiation, one should follow the graphs backwards, from right to 
left. Real non-zero minima of expansion are taken here to arise from the energy density in the 
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radiation field taking a negative value. Again, a cut-off in energy transfer may be needed to 
prevent this from occurring. 

5.3 A flat universe containing radiation and vacuum energy 

From (Tl3]) and (j!4p it can be seen that a spatially flat universe containing radiation (T = 4/3) 
and negative vacuum energy (7 = 0) evolves as 

a 2 c = c 2 Y (x) + c 3 vr J (x) (19) 

when s = s A, and as 

a 2 D = c 2 A, (y) + czBi (y) (20) 

when s = So, where x and y are the same as before. The form of these solutions are shown in 
figures [3] and HI 

Again, these solutions will take both positive and negative values for any set of c 2 , C3 and 
Ao, and likewise, solutions can be matched at the moment of crunch-to-bang to ensure that a 
remains real and positive. This is done in figures [3] and HI 




Figure 3: The evolution of the scale factor in a flat universe with radiation and vacuum ex- 
changing energy, as prescribed in equation ( HE) when < s oc A. 

Both of figures [3] and H] show universes which are initially oscillating, with scale factor 
oscillations of increasing amplitude that eventually end in an asymptotic period of continual 
expansion. This is due to the choice s > 0, which corresponds to energy being transferred from 
the radiation to the negative vacuum energy. Once again, the contrary case of energy transfer 
from vacuum energy to radiation can be considered by reversing the direction of time in these 
plots. 

In figure Owe have the energy density of the vacuum evolving as A = — |A |e~ 3 ' s °'*, so that 
as t — > 00 the energy density of the vacuum goes to zero. This has the effect of leaving a 
radiation-dominated universe (a 2 oc t) in the late-time limit, as can clearly be seen from the 
plot. The exchange of energy then becomes unimportant as subsequent evolution continues 
in the standard way, with the vacuum energy density having been forced to zero. At earlier 
times, however, the influence of the negative vacuum energy becomes increasingly important, 
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resulting in oscillations of ever- decreasing amplitude. In contrast to some of the previous cases, 
the energy density of the radiation is always positive in this scenario, because increasing the 
magnitude of the negative vacuum energy results in an increase in the positive energy density 
of the radiation. 




Figure 4: The evolution of the scale factor in a flat universe with radiation and vacuum ex- 
changing energy, as prescribed in equation ( \2(fy when < s =constant. 

In figure HI the transfer of energy from the vacuum to radiation occurs at a constant rate, 
resulting in an evolution of the vacuum energy density as A = — |A| + 3sot. For sq > 0, energy 
is transferred from the radiation to the vacuum, as in the previous example. Now, however, 
the energy density of the vacuum is allowed to become positive at sufficiently large t. This 
results in a period of rapidly accelerating expansion as t — > oo: not only is the evolution of 
this universe dominated by a positive cosmological constant, but the value of this "constant" 
is itself increasing with time. The early-time behavior is qualitatively similar in this example 
to the previous one: as t becomes increasingly negative, oscillations are increasingly damped. 
Once again, we are guaranteed a positive energy density of radiation, via the same mechanism 
as before. The principal difference here is that the decrease in amplitude of oscillations, as 
t decreases, is slower than before because the transfer of energy is now constant, and not 
proportional to A. 

Inclusion of non-zero spatial curvature in these models results in a more elaborate form 
of a(t), as given in equations f[T51) and f[T4"l) . However, although more difficult to express in 
a concise analytic form, the behavior of these solutions can be simply understood. In limits 
where radiation previously dominated the evolution (such as the t — > oo limit of figure [3]) spatial 
curvature now eventually dominates. Positive curvature causes a maximum of expansion, and a 
subsequent collapse, while negative curvature causes evolution towards a Milne universe, a act. 
In limits where the vacuum energy dominates expansion, the spatial curvature has a negligible 
effect, unless it is large enough to cause collapse before vacuum domination can occur. 

6 Discussion 

We have investigated the behavior of spatially homogeneous and isotropic oscillating universes 
when non-equilibrium behavior is present. We have found exact solutions for universes con- 
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taining multiple fluids and non-zero spatial curvature. These fluids include radiation and dust, 
scalar field, or vacuum energy. 

In the usual treatment, cosmologies with more than one fluid component are assumed to 
evolve with these fluids staying non-interacting. The energy of each fluid is then separately 
conserved, and the evolution of the scale factor can be straightforwardly obtained by solving the 
Friedmann equation. Here, we have considered a more general situation in which the different 
fluid components of the universe are allowed to interact. Total energy is still conserved in these 
models, but we have allowed energy to be transferred between the different components. This 
situation has often been discussed in the literature, with cosmological models in which entropy 
increases only at each moment of crunch-to-bang. Here, we have formulated the problem 
consistently and found exact solutions in which entropy production occurs realistically and we 
have used them to model a number of different possibilities. 

Firstly, we considered universes containing radiation and dust and universes containing ra- 
diation and scalar field. In these models the oscillatory nature of the solutions is due to positive 
spatial curvature. By allowing interactions between these components, we find oscillating so- 
lutions where the amplitude of successive cycles is allowed to vary. The case of radiation and 
dust is of obvious interest as these components are required for the nucleosynthesis of the light 
elements and formation of large-scale structure in the standard cold dark matter model. The 
case of radiation and scalar field is also of interest as scalar fields, if any exist in nature, will be 
influential in the early evolution of the universe. We find that the evolution of these universes is 
highly dependent upon the exact form of interaction between the fluids. For some interactions, 
such as those of the form considered in pQ, oscillating universe models are found to progress with 
each of their cycles having identical amplitude. For other interactions, it is found that energy 
transfer leads to cycles of varying amplitude. Again, the precise from of this variation is heavily 
dependent on the precise form of interaction. We find general methods for solving the Fried- 
mann equations, and investigate explicitly two particularly simple models. It has often been 
conjectured that oscillating anti-damped universes of growing amplitude may offer a solution 
to the flatness problem, without the need for inflation. By allowing the amplitude of each cycle 
to increase monotonically, it is suggested jH [5] that eventually the universe will end up being 
pushed closer and closer to spatial flatness, with each successive cycle being longer and longer 
lived. Our toy models show us that this may not be a particularly realistic expectation, even if 
there is no cosmological constant. We find that the increase in amplitude of the oscillations is 
halted in cases where the energy density of fluids is bounded from below: there is simply not 
enough energy to allow the cycles to become indefinitely large and long lived. In cases where 
the energy densities are not bounded from below, such as are often considered with scalar field 
models, we also find problems. When the energy density of these components become negative 
the universe experiences a non-zero minimum of expansion, instead of a crunch to singularity. 
This behavior is well known, and is often used to model non-singular oscillating cosmologies. 
However, we find that the continual transfer of energy required to create cycles of indefinite size 
and duration results in subsequent increases in the non-zero minimum of expansion. Thus, by 
sourcing energy from this negative energy scalar field it becomes more negative, and increases in 
the maximum of expansion are accompanied by increases in the minimum of expansion. This 
behavior does not appear to be consistent with a physically viable cosmology - that is, one 
which starts small and lives long enough for matter domination to occur. 
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In addition to these models, we also consider universes containing radiation and non-zero 
vacuum energy. The oscillatory nature of these solutions is either produced by positive spatial 
curvature, or by a negative vacuum energy. For spatial curvature to have any non-negligible 
effect on these cosmologies it must be strong enough to dominate the evolution of the universe 
before the effects of vacuum energy become significant, otherwise vacuum effects dominate. 
Again, we provide a general prescription for solving this problem for arbitrary functions pa- 
rameterising the exchange of energy. As before, we give two particularly simple toy models 
explicitly. In these simple models we neglect spatial curvature and allow a negative vacuum 
energy to collapse the universe. We find that if the vacuum energy is allowed to become positive 
then it quickly comes to dominate, and accelerating expansion rapidly ensues. On the other 
hand, if the vacuum energy becomes increasingly negative then the amplitude and duration of 
cycles is found to decrease as the vacuum energy comes to dominate, and collapses the universe 
increasingly rapidly. Finally, in examples where the exchange of energy causes the vacuum 
energy to approach zero, we find that the subsequent evolution of the universe progresses as 
a standard radiation-dominated Friedmann universe. Including the effects of spatial curvature 
complicates the precise analytic form of the solutions, but otherwise acts in the expected way. 
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Appendix: From s(a, t) to m(t 

It can be seen from (|2j) and ([3]) that we can obtain any invertible function m = m(r) by 
choosing the energy exchange parameter to take a form 

s=[a 3 y'(pa^)g(a)]- 1 , 

where / is the inverse function of m — m(t) and g(a)dT = dt defines the time coordinate r. 
Substituting this into d2J) gives 

riP^^p- = i, 

which integrates to 

f(pa^) = r-r 

or 

pa 31 = m{r - r ) = / -1 (r - r ), 

where r is a constant of integration, which can be trivially absorbed into a redefinition of the 
time coordinate r — > r + r . This gives us a prescription for the form of s(p,a) required to 
derive any invertible function m = m{r). 

We will often be interested in the simple forms 

mi = m e s ° T and m<i = m + s r (21) 
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where m and s are constants. It can now be seen from the above that these forms for m 
correspond to 

si = — and s 2 = (22) 

9 a 7 # 

respectively. 

These results can be used to find exact solutions for a(t) in universes with two interacting 
fluids, and non-zero spatial curvature. 
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